Normal heat conductivity in a strongly pinned chain of anharmonic oscillators. 
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We consider a chain of coupled and strongly pinned anharmonic oscillators subject to a non- 
equilibrium random forcing. Assuming that the stationary state is approximately Gaussian, we first 
derive a stationary Boltzmann equation. By localizing the involved resonances, we next invert the 
linearized collision operator and compute the heat conductivity. In particular, we show that the 
Gaussian approximation yields a finite conductivity k ~ J^W ^ or ^ ^he anharmonic coupling 
strength. 
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The analytical derivation of the Fourier's law of heat con- 
duction in extended Hamiltonian systems remains to this 
day a challenging question. See 0, for a review of 
various aspects of the problem. A qualitative explana- 
tion for Fourier law has been provided by Peierls for the 
case of crystalline solids. Peierls idea essentially con- 
sists in considering the lattice excitations of the crystal, 
the phonons, as particles interacting weakly through the 
nonlinear forces acting between the atoms. By analogy 
with the case of dilute gases, a Boltzmann equation for 
honons has been derived and refined over the years, see 
and references therein. This approach yields a qual- 
itative comprehension of the heat transfer in terms of 
collisions between phonons. The explicit computation of 
the heat conductivity amounts next to invert the colli- 
sion operator between phonons. To this day, this com- 
putation has not been performed. This is due to the fact 
that the localization of the resonances in the combination 
of frequencies of the phonons involved in the collision is 
a notoriously delicate problem. In this Letter, we solve 
this problem in the case of a strong harmonic pinning. 

Crystalline solids submitted to a heat flow are often 
modelized by a chain of oscillators where each oscillator 
move around some equilibrium position located on a reg- 
ular lattice. Each oscillator interact with its nearest neig- 
bours and also possibly with some substrate. The chain is 
coupled to heat baths at different temperatures by some 
random forcing and friction acting only on the particles 
at the boundaries of the chain. Although mathematical 
properties of such systems have been extensively stud- 
ied, see 6] for a review, very little is known analytically 
about their physical properties, except in cases where 
the interactions are quadratic, pj. Numerous numerical 
simulations have been performed, see p| and references 
therein, strongly supporting that Fourier law is verified 
in pinned anharmonic chains. In particular, the conduc- 
tivity of our model has been studied numerically Q, 3 
and shown to be finite, although with a different scaling 
in temperature than what the Gaussian assumption of 
the Peierls theory predict, as we will see below. 

As we want to concentrate on the analysis of the oper- 
ator of collision between phonons, we derive a stationary 
Peierls-Boltzmann equation for a translation-invariant 



chain. The equation for the evolution of the average en- 
ergy current </& carried by the mode of wavenumber k of 
the chain will have the form, 



dt 



J kit) = ~lJk[t) + AA4(«) + S k 



(1) 



The first term on the RHS of this equation is a damp- 
ing term coming from the coupling of the chain to an 
external friction represented by the coefficient 7. The 
second term comes from the Hamiltonian evolution and 
vanishes when the anharmonic part of the interaction is 
zero, i.e., when A = 0. The third term represents an 
external "creation" of current at constant rate. It mod- 
elizes the local temperature gradient in the non-periodic 
chain. Equation will therefore yield a balance equa- 
tion for the current carried by the mode k in the sta- 
tionary state. We will show that under the assumption 
that the stationary state is approximately Gaussian, J\fk 
is analogous to a friction term, yielding in the end a con- 
ductivity k ~ X i Ti . From a mathematical point of view, 
there is an important difference between the part coming 
from the external friction and the part coming from A4 . 
Whereas the friction acts as a diagonal operator on the 
vector J, we will see that A4 ~ A(J)k where A is a non- 
diagonal linear operator. This is the collision operator 
of the linearized Boltzmann equation for phonons. The 
properties of this operator are essentially governed by res- 
onances in the four-body interactions between phonons, 
which we analyze in the limit of a large harmonic pin- 
ning strength fi. In particular, we show explicitly that 
the origin of the finite conductivity is to be found in 
/^-dependent resonances corresponding to umklapp pro- 
cesses. As we concentrate only on the effect of the term 
coming from the nonlinearity A4, the source term Sk 
is put by hand. However, in the non-translation invari- 
ant case, it comes naturally from the transport term of 
the equations and is roughly given by uj~ 1 (k)\7uj(k)\7T 
where u(k) is the dispersion relation and VT is the local 
temperature gradient, see 

The Model. We consider particles located on a periodic 
lattice described by their positions and momenta coordi- 
nates (qi,Pi) with i G 1, . . . ,N. For the sake of notational 
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simplicity, we assume that N is divisible by 4. The dy- 
namics of the particles is governed by the Hamilton equa- 
tions and an external random forcing. The Hamiltonian 
that we consider is of the form 



N 
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The anharmonic interaction is only in the on-site part 
of the Hamiltonian. The same analysis as found below 
may be carried out when the interaction between nearest 
neighbours contains also a quartic part. The fact that 
there is an on-site (large) harmonic interaction makes 
the distinction between the two cases irrelevant. We first 
consider the harmonic case A = 0. We introduce the 
Fourier coordinates for the periodic linear chain by Qk = 
737 EjLi e^ ^with -N/2 + 1 < k < N/2. The P k 
coordinates are defined in a similar fashion. We recall for 



further purposes that Q k 
complex coordinates, 

At- 



-k and Pi 



(P k ± iuJ k Qk), 



P- k - In the 



(3) 



with lo\ 



u. > (fi +4 sin (^-)), the Hamilton equations 
for the linear periodic chain read 



dAf = ±iuj k Afdt. 



(4) 



Those equations give the temporal evolution of the ampli- 
tudes of the waves with wave number k traveling through 
the chain in the positive or negative direction. Our (non- 
equilibrium) model is defined in the linear case by the 
system of stochastic equations 

dAf = ±iu> k Afdt -\{A+ + A% )dt + dW£ , (5) 
where the Wiener processes W k satisfy the relations 



7T 



+ (s-s')Ta(^r) <Si/,-:-A-')'//. 



(6) 

In (JSJ) , the term in 7 is the usual friction term acting only 
on the momenta Pk . The Wiener processes represent the 
random injection of energy in the system. The first term 
in © is the usual term giving the equilibrium dynam- 
ics. In the second term of (JBJl, a is an odd, dimensionless 
function of period it. Ultimately, the choice of physically 
relevant a is fixed by the effect of the local temperature 
gradient in the transport term of the non-translation in- 
variant chain. Finally, we interpret r in the second term 
of JHJ as the non-equilibrium parameter of our model. 
When r is zero, the periodic chain is in thermal equilib- 
rium at temperature T. When r^O, a stationary heat 
current proportional to r is created in the chain. This 
is achieved, via the factor s — s' , by exciting more the 
waves travelling in one direction than the ones travelling 
in the opposite direction. 



In order to see those points, we recall first that the 
energy current in the periodic chain reads 

J=— Yl sin(2^fc/iV)3(Q_ fe P fe ), (7) 



fe=-JV/2+l 



where Q(z) denotes the imaginary part of a complex z. 
We define the two-point correlation functions in the sta- 
tionary state, 



(8) 



where (.) denotes the average w.r.t. the stationary 
measure. Since 3((Q_ fc P fe )) = |($£l fe - ®l + - k ) = 
i«|A+| 2 )-(LA-| 2 )),weget 



2 N l 2 

fc=-AT/2+l 



in(2nk/N) [<| A+\ 2 ) - (\A^ | 2 )] . 

(9) 



This identity expresses the fact that the energy current is 
the difference of the current carried by waves traveling in 
the positive and negative directions. From (jSJ, we derive 
the equation for the two-point correlation functions, 

i(su> k + «w)*$ = + ~ ^ + \iki 

(10) 

with If'ti given in (JSJ. In the case where r = 0, (|10|l is 
readily solved to yield, using ©, 



k,k' 



— S(k + k')S(s + s'). 
^k 



(11) 



This is of course equivalent to (|P fc | 2 ) + u 2 k {\Qk\ 2 ) = 2T 
and (P k Qi) =0,Vk,l. 

When t ^ and for I = — k, one obtains from (|10|l and 
©, 

Jk = *fck-*i;+k = *n- 1 <*C£)> ( 12 ) 

which is the balance equation © in the stationary state 
for the harmonic case. We note that by definition, 
Jfc = — J-fc. Using © and the fact that a is an odd 
function, we get (J) = 2lq 2 ^~ 1 t J^ 2 a(x) sm(2x)dx for 
large N. As we interpret r as a local non-equilibrium 
parameter, the conductivity in the harmonic case, n = 
(J) /t ~ 7 _1 w 2 , behaves as in the harmonic chain cou- 
pled to self-consistent reservoirs 2J. It reflects the fact 
that when A = in the balance equation QJ, the only 
term damping the current carried by each mode comes 
from the external friction. We will see below that the an- 
harmonic interactions modify that behaviour of the con- 
ductivity. 

Let us now consider the anharmonic case A ^ 0. The 
equations of motion JSJ become 



dAt 



±tu k A±-l(Aj 



~ A k)- 



iX 
N 



dt+dW*, (13) 



3 



where 

Rk = y y <5(fc - fci - /s 2 - fc 3 ) sis 2 s 3 

ki,k 2 ,k 3 si,S2,S3 

w^M 2 ^, (14) 

and L kklk2k3 = (16uj k u; kl uj k2 uj k:s y 1/2 . Here and below, 
unless otherwise specified, the sums are over the ki such 
that -JV/2 + 1 < ki < N/2 and s, = ±1. 

From (|13fl . the equations for the n-point correlation 
functions in the stationary state, 



(15) 



read 



jo) (")$(") = 7r (n) ($ (n) )+i— -M (n) ($ ( " +2) )--!-/ (n) ($ ( ™ _2) ), 

N ' 2 

(16) 

where is the combination of frequencies, 

n 



and 



(r( n) ($<»>))£;;;£ =£(*<»>)£;;;;: 



^E( $(n) )" 1 ;::X'- ,Sn - ( 17 ) 



I*™- 1 gathers the effects of the random forcing, the special 
case 1^ being given in ©. The explicit expression of 
MW will be given below in the relevant cases. For n = 2 
and fci = — k 2 = fc, one gets from (II 61) the balance equa- 
tion in the stationary state, proceeding as to obtain (|12H , 

7*4 + ^ = 4™^), (18) 

where J k has been defined in (|12Jl and 

*^(Af( 2 )(<i>( 4 )))+- fe -(M( 2 )(<i>( 4 )))-+ fe 

2i > > S(k + ki + k 2 + k 3 ) sis 2 s 3 



ki,k 2 ,k 3 si,S2,S3 



with 



The remainder of this paper is devoted to solving equa- 
tion (ISJ for J fc . 

Closure and Linearization. In this section, we derive 
an explicit expression for 'J in terms of J. The first step is 
to approximate the n-point correlation functions l|15ll by- 
products of <I>( 2 ), turning (fl%)) into a close but nonlinear 
equation. This Gaussian closure assumption on the non- 
equilibrium stationary measure is widely used in various 



forms, such as in, e.g., the study of dilute gases. The 
second step is to linearize the resulting equation around 
the equilibrium solution ifTTJl . This is done in the spirit 
of taking a small temperature gradient or, in our case, a 
small parameter r. 

Closing as described above the expression l|19|) for ^ 
yields zero after linearization, so that equation 1)18(1 is 
identical to equation i|12|) for the harmonic case. We 
thus first express <&( 4 - ) in (|19|l in terms of <f>( 6 ). Inverting 
(PI for $( 4 ), one can show that at lowest order in 7 



c >[/'*( 4+ ) > i +SlS2S3 
)k,kiM,) 

where 

tri \sis2s3s4 

^ Si 7/ki,ka,fc3,foi — 



A 



S2S3 

ki,k 3 



(21) 



20 



47 2 + (s^i + s 2 uj k2 + s 3 uj k3 + s 4 w fe4 ) 2 



In l(2ip. we have dropped a term containing 3[$( 6+ )] as 
well as the contribution from the Ito term since 
they vanish later in the computation. We note that when 
7^0, 

( n i)l\ S MMM -> <^( s i^i +s 2 uJ k2 +s 3 uj k3 +s A uj ki ). (22) 

We next plug <|2"T|l into ifT^I and use the Wick formula 
to express $( 6 - ) in terms of Linearizing around the 

equilibrium solution (|ll|l results in an expression for W 
which depends solely on J and contains terms which may 
be represented graphically by 12 topologically distinct 
diagrams. Those diagrams are structurally similar to the 
usual diagrams of the equilibrium A(/> 4 theory. The mam 
difference comes from the presence of the quantity which 
is the unknown of the problem, i.e., the vector J. 

For large pinning interaction fj,, the result of the com- 
putation is (here and below, c denotes combinatorial pos- 
itive factors), 



cT 2 iX^,^ 

' l.n 



+ H 

l>k,l,n,k+l+ 



n [J k + Jl+ J n — Jk+l+n] ■ 

(23) 

The RHS is the linear operator A mentioned in the in- 
troduction, that is, the linearized collision operator of 
the Boltzmann equation. We remark that, in quantum- 
mechanical language, only the interactions conserving 
the number of phonons contribute, as the other inter- 
actions are free from resonances. We will estimate l|23[) 
in the limit of large N and for small 7, more precisely, 
1 ^> 7/cj 3> 1/N. We first observe that the main contri- 
butions to the sum will arise from the resonances in fly, 
i.e., for lattice points (I, n) near the zeros of the function 

Jx ■ [ 2 ' 2 J 



My, a) = \/m 2 + 4 sin 2 (x) + n 2 + 4 sin 2 (y) 

- J u 2 + 4sin 2 (z) - \J /i 2 + 4sin 2 (a; + y + z). 



(24) 
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A careful analysis reveals that for /i 2 > 0, the zeros of 
f x form three smooth curves. Two of these curves are 
obvious and given by x + z = 0, and y + z = 0. These 
resonances do not contribute to the sum, however, since 
the combination of J's in (|23[l vanishes for k + n = 
or I + n = 0. This corresponds to the so-called normal 
processes. The third curve, corresponding to the umk- 
lapp processes, depends on fi and is difficult to localize 
explicitly. For large /i, it is given by x + y — §+0(1/ u 2 ). 
Performing the sum over I in l!2.'iH in the above-mentioned 
parameter regimes thus yields at lowest order in 7 and 
1//J (resonances are at I = — k + N/2), 



to 7 /I 5 ^ 



sm (jv ) cos ( 



7r(k— n) 1 

N • 



(25) 



where we have defined JjT = \{Jk — Jk+N/2)- Before 
proceeding to solving equation (|18f) . we note that any 
vector J can be decomposed as J — J + + J~ where J~ 
\{Jk + Jk+N/2) have the symmetry 



as above and j£ 



properties J^ +N / 2 = ^^k' ^ then follows from 125fl that 
J + does not contribute to "J, and that ^ = "f^. 

Solution to the Current Equation. Recall that J 
is odd and periodic of period N. We first observe that 
J + does not contribute to the average current (J) given 
by ©. This follows from £ fe sin(27rfc/A) J+ = 0. Since, 
in addition, J is mapped by (|25|l into a vector \& with 
^ = and since J + does not contribute to 'J, we need 
only to consider in equation (|18|1 odd forcings a and odd 
currents J satisfying a(^ + f ) = — a(^) and J = J~ . 
We denote by S~ the subspace of such vectors J, in which 
equation l(TH|) becomes 



jJ k + c ^5 C{J) k = *raA , 
oj' ii J N 



where C : S — ► <S is given by 



jfe + jfe 



^^|sin(^))cos(^); 



(26) 



(27) 



We now proceed to analyze the linear operator C. The 
subspace S~ has dimension A/4, and a basis for S~ is 



given by J/} = sin(27r(2n+l)|0, n = 0, . . . ,N/4- 1. We 
let ^4" = C(J n ). A direct computation shows that the 
set of ^4™, n = 0, . . . , N/4 — 1, also forms a basis of S~ 
and that £ is uniformly invertible (in N). This implies 
that the first term on the LHS of (|26|l is negligible for 7 
small. Furthermore, it follows from sva.{^j^-)J^ = 
for n ^ that only contribute to the current 
Therefore, the only contribution of the noise a to the 
current is the component of a along A , say a , where 



Al = £(J°)k = 2sign(fc) - - sign(fc) 



N 



(28) 



One thus finally obtains from ©, 



{J) 



1 i9,,5 

■ a t , 



(29) 



A 2 r 2 

Conclusions. As mentioned in the introduction, the 
result of @ shows that the local gradient of tempera- 
ture imposes the choice a(k) ~ a; (fc)Vw(fc), that is, 
/i~ 2 sin(27rfc/iV) for large /i. It has a non-zero compo- 
nent along A and the corresponding conductivity is thus 
given by k = (J) /t ~ jsj^- O ur explicit treatment of 
the resonances and the inversion of the linearized Boltz- 
mann operator allows to justify rigorously the physical 
picture that the mechanism responsible for the normal 
conductivity in pinned anharmonic chains is the four- 
body umklapp collisions between phonons. Indeed, con- 
sidering a cubic instead of a quartic interaction in the 
Hamiltonian would yield an expression analogous to i|22|) 
but with a combination of only three frequencies. The 
latter is always non-zero when /1 =/= 0, leading to an in- 
finite conductivity. We finally point out that when the 
pinning goes to zero, the localization of the resonances 
and the inversion of the linearized collision operator are 
much more difficult and there is no reason for expression 
(EHJ to be valid. 
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